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ABSTRACT 

HCN luminosity is a tracer of dense molecular gas, n(i?2)~3 x 10 4 cm -3 , associated with star- 
forming giant molecular cloud (CMC) cores. We present the results and analysis of our survey of 
HCN emission from 65 infrared galaxies, including nine ultraluminous infrared galaxies (ULIGs, 
Lir^1O 12 L ), 22 luminous infrared galaxies (LIGs, lO n L < Lir^10 12 L q ), and 34 normal 
spiral galaxies with lower IR luminosity (most are large spiral galaxies). We have measured the 
global HCN line luminosity, and the observations are reported in Paper I. This paper analyzes 
the relationships between the total far-IR luminosity (a tracer of the star formation rate), 
the global HCN line luminosity (a measure of the total dense molecular gas content), and the 
CO luminosity (a measure of the total molecular content). We find a tight linear correlation 
between the IR and HCN luminosities Lir and Lhcn (hi the log-log plot) with a correlation 
coefficient R = 0.94, and an almost constant average ratio Lir/Lhcn = 900L Q /Kkms _1 pc 2 . 
The IR-HCN linear correlation is valid over 3 orders of magnitude including ULIGs, the most 
luminous objects in the local universe. The direct consequence of the linear IR-HCN correlation 
is that the star formation law in terms of dense molecular gas content has a power law index of 
1.0. The global star formation rate is linearly proportional to the mass of dense molecular gas 
in normal spiral galaxies, LIGs, and ULIGs. This is strong evidence in favor of star formation 
as the power source in ultraluminous galaxies since the star formation in these galaxies appears 
to be normal and expected given their high mass of dense star-forming molecular gas. 

The HCN-CO correlation is also much tighter than the IR-CO correlation. We suggest 
that the nonlinear correlation between Lir and Leo may be a consequence of the stronger 
and perhaps more physical correlations between Lir, and Lhcn and between Lhcn and Leo- 
Thus, the star formation rate indicated by Lir depends on the amount of dense molecular gas 
traced by HCN emission, not the total molecular gas traced by CO emission. One of the main 
arguments in favor of an active galactic nucleus (AGN) as the power source in ULIGs is the 
anomalously high ratio Lir/Lco or Lir/M(H2) or high star formation rate per M Q of gas, 
compared with that from normal spiral galaxies. This has been interpreted as indicating that 
a dust-enshrouded AGN is required to produce the very high luminosity. Viewed in terms of 
the dense gas mass the situation is completely different. The ratio Lir/Lhcn or LiR/Mdensei 
a measure of the star formation rate per solar mass of dense gas is essentially the same in all 
galaxies including ULIGs. The ratio Lm/Mdense * s virtually independent of galaxy luminosity 
and on average LiR/Mdenso ~ 90L Q /M Q , about the same as in GMC cores but much higher 
than in GMCs. We find that ULIGs simply have a large quantity of dense molecular gas and 
thus produce a prodigious starburst that heats the dust, produces the IR, and blocks all or most 
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optical radiation. The HCN global luminosity may be used as an indicator of the star formation 
rate in high-redshift objects including hypcrluminous galaxies. 

The HCN/CO ratio is an indicator of the dense molecular gas fraction, and gauges the 
globally averaged molecular gas density. We find that the HCN/CO ratio is a powerful 
starburst indicator. All galaxies in our sample with a high dense gas mass fraction indicated 
by Lhcn/Lqo > 0.06 are LIGs or ULIGs. Normal spirals all have similar and low dense gas 
fractions Lhcn/^co = 0.02 to 0.05. The global star formation efficiency depends on the fraction 
of the molecular gas in a dense phase. 

Subject headings: galaxies: ISM — galaxies: starburst — infrared: galaxies — ISM: molecules 
— radio lines: galaxies 

1. INTRODUCTION 

Stars are born in the molecular interstellar medium, the raw material for star formation. In the Milky 
Way, all star formation essentially takes place in molecular clouds and most star formation takes place 
in giant molecular clouds (GMCs; Solomon, Sanders, & Scoville 1979) with mass M > 10 5 M Q and not 
the diffuse neutral ISM dominated by atomic hydrogen (extended HI gas disk). The star formation rate 
(SFR) of molecular clouds can be estimated from the far-infrared (FIR) luminosity emitted by the warm 
dust heated by embedded high-mass OB stars (e.g., Mooney & Solomon 1988). The mass of molecular 
gas can be determined from the CO luminosity calibrated by 7 ray flux from the interaction of cosmic 
rays with hydrogen molecules (e.g., Bloemen et al. 1986) or by dynamical cloud masses determined from 
CO kinematics for virialized individual molecular clouds (Solomon et al. 1987; Young & Scoville 1991). 
These methods are in good agreement (see Solomon & Barrett 1991). All strong high-mass star formation 
regions are associated with GMCs, especially the cores of GMCs. The ratio of FIR luminosity to the CO 
luminosity or to the cloud mass, a measure of the SFR per solar mass of the cloud and an indicator of star 
formation efficiency (SFE), ranges over a factor of 100 for different clouds, and over a factor of 1000 from 
clouds to the cores of GMCs (e.g., Mooney & Solomon 1988; Plume et al. 1997). An understanding of the 
physical conditions in GMCs and their relation to galactic dynamics is a prerequisite to the understanding 
of the star formation process, the SFR in galaxies and starbursts. 

Star formation in galaxies is closely tied up with the local gas density, as formulated in the Schmidt 
(1959) law, although the important component is the molecular gas. Globally, the SFR correlates with the 
molecular gas content in galaxies, as traced by CO emission, including luminous and ultraluminous infrared 
galaxies (LIGs and ULIGs 1 ). In Galactic star- forming regions, active high-mass star formation is intimately 
related to the very dense molecular gas in the cores. While the canonical molecular gas tracer CO shows 
strong emission in cloud cores, it is not specific enough to reveal their star formation potential. The bulk 
of the cloud material traced by CO observations is in the GMC envelopes and is at a much lower density. 
The physical conditions of active star-forming GMC cores are better revealed by emission from very high 
transition CO lines (in the submillimcter regime) and high dipole-momcnt molecules like CS and HCN. 



1 LIGs: 10 12 L Q >L IR > lO n L , ULIGs: L IR >1O 12 L (to be exact, 1O 119 L in this paper). For a definition of the total IR 
(8 to 1000/^m) luminosity Ljr and FIR luminosity I/fir, see Sanders & Mirabel (1996). The value of Lir is generally larger, 
by up to ~20% as it includes both 12 and 25/j.m emission, than Lpir. However, we often simply refer the total IR emission as 
FIR in this paper. 
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HCN is one of the most abundant high dipole-moment molecules that traces molecular gas at densities 
n(£/2)~3 x 10 cm" 3 , more than 2 orders of magnitude higher than that traced by CO (^,300cm -3 ). Many 
HCN (1-0) observations have already been conducted by different groups (Nguyen-Q-Rieu et al. 1989, 1992; 
Henkel et al. 1990; Solomon, Downes, & Radford 1992; Israel 1992; Heifer & Blitz 1993; Aalto ct al. 1995; 
Curran, Aalto, & Booth 2000) in a variety of nearby galaxies. These previous observations contain small 
samples with frequent overlap in the sample selection. The total number of galaxies detected in HCN is still 
small (,$30), and for many galaxies only the central position was observed. Gao & Solomon (2003, hereafter 
Paper I), in the companion paper, have presented a systematic survey of global HCN luminosity that more 
than doubles the number of galaxies (~ 60) observed in HCN. 

Although the dense molecular gas is strongly concentrated in the central regions (,$,lkpc), HCN 
mapping of a dozen nearby galaxies out to a diameter of ~ D25/4 (Gao 1996, 1997; Y. Gao & P.M. Solomon 
2004, in preparation) shows that a substantial fraction of HCN emission originates from the inner disks 
outside the central <~ 1 kpc. All previous HCN observations in external galaxies including a few HCN 
maps (e.g., Nguyen-Q-Rieu et al. 1992; Reynaud & Downes 1997; Heifer & Blitz 1997a) have primarily 
been observations of the galactic central regions. In most cases, the total HCN emission from the entire 
galaxy has not been measured. The recent HCN observations in 20 Seyfert galaxies are primarily of the 
galaxy centers, and many are nondetections or marginal detections (Curran et al. 2000). Confusion may 
occur when the results drawn from the HCN observations of central regions of galaxies are compared with 
the global properties of galaxies. Aalto et al. (1995) have not found a correlation between the molecular 
line intensity ratio, /co/^hcn, and FIR emission, or measures of star- forming activity in their sample of 
10 interacting galaxies, which appeared to be in conflict with the original findings of a tight FIR-HCN 
correlation in another sample of 10 LIGs/ULIGs and spiral galaxies where the total HCN emission was 
measured (Solomon et al. 1992). The situation has now been clarified by our new HCN survey of ~ 60 
galaxies (Paper I), which measured the global HCN emission in a wide range of galaxies. There is indeed a 
tight FIR-HCN correlation in a statistically significant HCN sample. 

It is clear from all previous observations that the molecular gas in the central regions of spiral galaxies, 
starbursts, and LIGs/ULIGs is much denser than the molecular ISM in the disks of spiral galaxies (e.g., 
Nguyen-Q-Rieu et al. 1992; Solomon et al. 1992; Heifer & Blitz 1997a; Wild & Eckart 2000). From various 
observational studies of star-forming regions of GMCs in the Milky Way, it is also clearly shown that the 
active star-forming regions in the disk are the dense molecular cloud cores (e.g., Mooney & Solomon 1988; 
Plume et al. 1997; Pirogov 1999; Evans 1999) rather than the entire molecular envelopes of GMCs. The 
SFE of active star-forming clouds that are associated with IR sources readily apparent on IRAS 60 and 
100 fim images can be 100 times higher than the IR-quiet clouds with no apparent IR sources revealed by 
IRAS (Mooney & Solomon 1988). Star formation efficiency (Lir/£co) can var y over a factor of 100 as well 
in galaxies. 

Normal spiral galaxies have an SFE similar to that of Galactic GMCs, whereas the SFE of ULIGs/LIGs 
and advanced mergers can be more than an order of magnitude larger (e.g., Solomon & Sage 1988; Solomon 
et al. 1997; Gao & Solomon 1999). These differences in SFE can be understood in terms of the different 
dense molecular gas content as traced by HCN observations. 

For 10 galaxies, including both LIGs/ULIGs and normal spiral galaxies, Solomon et al. (1992) show 
that there is a tight correlation between the ratio of the IR and CO luminosities £ir/Lco and the HCN/CO 
luminosity ratio I/hcn/-^co in addition to the excellent correlation between Lir and Lhcn- This is now 
fully confirmed with our current HCN study of a sample of 65 galaxies (§3). These correlations demonstrate 
a close relationship between the SFR and the dense molecular gas reservoir in galaxies. The SFE depends 
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on the fraction of available molecular gas in a dense phase (Lhcn/£co)> and the dense molecular content 
of even gas- rich spirals is much less than that of LIGs/ULIGs of comparable total molecular gas content 
(Solomon et al. 1992; Radford 1994). Since CO emission traces most of the molecular gas mass and is not 
necessarily a specific tracer of the dense molecular gas (e.g., Mauersberger & Henkel 1993; Evans 1999) or 
the IR luminosity from star formation (Mooney & Solomon 1988), CO alone can give a misleading picture 
of the densest molecular gas in a galaxy. 

In Galactic plane GMCs, for example, essentially all OB star formation occurs in the cores of GMCs 
with strong CS and HCN emission. The ratios of CO/CS and CO/HCN intensities for Galactic disk GMCs 
are much larger than for Galactic center clouds (Lee, Snell, & Dickman 1990; Jackson et al. 1996; Plume 
ct al. 1997; Heifer & Blitz 1997b) and an order of magnitude larger than for the archetypal ULIG Arp 220 
(Solomon, Radford, & Downes 1990). In some galaxies, the molecular gas in the center is much more 
prominent in HCN emission than in CO. A good example is the center of the Seyfert 2/starburst hybrid 
galaxy NGC 1068, where all interferometric maps demonstrate that the nuclear region is more prominent 
than the rings or spiral arms when viewed in HCN, while the opposite is true in CO emission (Tacconi 
et al. 1997, 1994; Heifer & Blitz 1995; Jackson et al. 1993). A similar trend is also observed in the centers 
of M51 and NGC 1097 (e.g., Kohno et al. 1996, 2003). 

Here we utilize a large, statistically significant sample with observations of global HCN emission from 
65 spiral galaxies, LIGs, and ULIGs — 53 from the systematic HCN survey of Paper I, 10 from Solomon, 
Downes, & Radford (1992), plus two from the literature. These galaxies range over 3 orders of magnitude 
in FIR luminosity. We analyze the various relationships among the global HCN, CO, and FIR luminosities. 
We further discuss the physical relationship between the dense molecular gas content and the rate of 
high-mass star formation in galaxies. 

The HCN sample and observations are briefly reviewed in §2. Section 3 presents the results. Section 
3.1 is a comparison of the IR-HCN and IR-CO correlations. Section 3.2 concentrates on the importance 
of the dense gas mass fraction and the star formation efficiency. Section 3.3 and the Appendices present 
multi-parameter fits to the data and the effect of dust temperature. Section 3.4 has a brief summary of all 
results. In §§4.1-4.3, we discuss the importance of HCN as a tracer of star-forming molecular gas, the star 
formation rate as a function of the dense gas mass, and the global star formation law. Section 4.4 discusses 
the HCN/CO ratio as a starburst indicator. Section 4.5 discusses the origin of FIR emission from spirals 
and ULIGs. Section 4.6 briefly speculates on the implications for hypcrluminous infrared galaxies at high 
z. Finally, we summarize the main points of our study in §5. 

2. THE HCN SURVEY 

The detailed descriptions of our survey sample and HCN (and CO) observations were given in the 
companion paper (Paper I). Our HCN survey sample is drawn from samples in recent CO observations of 
galaxies showing strong CO emission (e.g., the CO antenna brightness temperature much larger than 100 
mK for normal spiral galaxies and larger than 20 mK for LIGs/ULIGs). All truly IR- bright galaxies with 60 
or 100/zm emission larger than 50 or 100 Jy, respectively, have also been included. Essentially all galaxies 
with strong CO and IR emission in northern sky have been chosen for the HCN survey. 

We carried out several observing runs mainly with the former NRAO 2 12m telescope at Kitt Peak, 



2 The National Radio Astronomy Observatory is operated by Associated Universities, Inc., under cooperative agreement with 
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the IRAM 30m telescope at Pico Veleta near Granada, Spain, and the FCRAO 14m telescope for most of 
our observations. We here only mention our observing strategy of using different telescopes. The IRAM 
30m was mostly used to observe ULIGs (rather distant) and to map out some nearby starburst galaxies 
of smaller optical diameters, given the small matching beam size as compared with the source extent for 
effective mapping. Essentially all other observations were conducted with the NRAO 12m so that one beam 
measurement can cover almost all HCN emission from relatively distant galaxies, including some merging 
galaxy pairs. Also almost all nearby large spiral galaxies have been at least mapped along the major axes 
with the NRAO 12m, and only a few of them were initially tried mapping with the QUARRY receivers at 
the FCRAO 14m (because of its limited sensitivity). 

The goal of the HCN survey of Paper I is to determine the total HCN emission from the whole inner 
disks or the entire galaxies in a large sample of galaxies with a wide range of IR luminosity. Combined with 
the HCN data of 12 galaxies, 10 from Solomon et al. (1992), including a half-dozen LIGs/ULIGs, and the 
other two galaxies, M51 (Nguyen- Q-Rieu et al. 1992) and NGC 4945 (Hcnkel, Whiteoak, & Mauersberger 
1994) that were mapped extensively in HCN, we have a large, statistically significant sample of 65 galaxies 
with globally measured HCN emission. In addition, there are more than a dozen nearby large spiral galaxies 
with HCN detections toward the galactic nuclei (Nguyen-Q-Rieu et al. 1992; Heifer & Blitz 1993; Aalto 
et al. 1995; Curran et al. 2000), including LIGs (e.g., NGC 3256, Casoli, Dupraz & Combes 1992) that do 
not overlap with our HCN survey sample 3 . In principle, we could further enlarge our HCN sample to a 
total of about 80 galaxies. These observations of nearby galaxies did not measure the total HCN emission, 
however, with only the central/nuclear HCN emission observed. Therefore, most HCN data in these nearby 
galaxies can thus only be used to set lower limits to the total HCN luminosities. Further HCN observations 
at least along the major axes are required to map out the total HCN emission. For the sake of completeness, 
consistency, and uniformity of our sample, these galaxies with lower HCN limits in the literature are not 
included in our sample. 

The derived global properties of the line luminosities and various luminosity ratios of galaxies in the 
HCN survey sample are listed in Table 1 (cf. Table 4 in Paper I), together with a dozen galaxies where the 
total HCN luminosities are available or can be estimated from the literature. There are also several other 
fairly distant galaxies in Curran et al. (2000) and Aalto et al. (1995) in which the global HCN luminosities 
were supposedly measured, but we did not include them in our sample as their data seem to have low 
signal-to-noise ratio. Mrk 273, one of the galaxies that overlaps in both of our sample and that of Curran 
et al. (2000), for example, was claimed to be detected with an extremely large HCN/CO luminosity ratio 
of 1. Our high-quality IRAM 30m spectra (Paper I), however, clearly indicate an HCN/CO intensity ratio 
of ~ 1/6 and a HCN/CO luminosity ratio of 0.23, which is consistent with all other ULIGs that have been 
significantly detected in HCN so far, rather than an unrealistic ratio of <~ 1. 

We also have some limits to the total HCN line emission owing to insufficient mapping of HCN (about 
10 % of the galaxies in our sample), but our HCN lower limits are probably close to the true values as we 
have some off-nucleus HCN measurements besides the central beams on the nuclei in most cases. We keep 
these sensitive HCN limits in Table 1 and include these data for the various analyses. 



the National Science Foundation. 

3 However, Curran et al. (2001) submitted and published several HCN maps (though mostly still very limited spatial coverage 
of only central l'-2' regions) of nearby Seyfert galaxies after we submitted our paper. 
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3. RESULTS AND ANALYSIS 

3.1. Comparison of the IR-HCN and IR CO Correlations 

The principal observational result from this survey is the tight linear relation observed between 
far-infrared luminosity, Lir, and HCN line luminosity, Lrcn, shown in Figure la. This very good 
correlation extends over 3 orders of magnitude in luminosity and includes normal spiral galaxies and 
luminous and ultraluminous infrared galaxies (LIGs/ULIGs). Figure lb shows the correlation between Lir 
and Leo of the same sample, which has a larger scatter than the FIR HCN relation and most importantly 
steepens with higher IR luminosity, showing the well-known result that LIGs/ULIGs, although rich in 
molecular gas, have a substantially higher IR luminosity per unit CO luminosity or per solar mass of 
molecular gas H 2 (e.g., Solomon et al. 1997). 

A least-squares fit using all the data but excluding HCN limits yields a power- law slope of 1.00 ± 
0.05 and 1.25 ± 0.08 for Lir-Lhcn and Lir-Lco, respectively. The corresponding correlation coefficients 
(squared) are R 2 = 0.88 and 0.77. The slopes change to 1.05 ± 0.05 and 1.44 ± 0.08 for Lir-Lhcn and 
Lm-Lco, respectively, if an orthogonal fit is used. Including the sources with HCN limits in the fit has 
little effect on the slope or the fit. The best-fit (logarithmic) relation between HCN and IR luminosities 
excluding galaxies with upper or lower limits is 

logL m = 1.00(±0.05)Zo.gL HC N + 2.9, or Lir/Lhcn = gOOLo/Kkms^pc 2 , (1) 

indeed a linear relation. When we use all galaxies including HCN limits the correlation remains almost the 
same with logLm = 0.97^ogLucN + 3.1 and with the same correlation coefficient (R=0.94). 

At first glance of Figure 1 , the correlation between Leo and Lir may appear to be nearly as good as 
that of the Lrcn-Lir correlation, even though there is significant difference in the correlation coefficients 
and the dispersion from the fit. The most obvious difference occurs at the high Lir end. In Figure lb, a 
line of constant Lir/Lco = 33L Q /Kkms _1 pc 2 fits the lower luminosity galaxies but lies well below almost 
all (28/31) of the high-luminosity galaxies (LIGs/ULIGs). The correlation between Lir and Lhcn still fits 
the ratio determined from lower IR luminosities (Lir < 1O 11 L ), but the correlation between Leo and 
Lir does not. (The lines shown in Fig. 1 are the fits of a fixed slope of 1 to low-luminosity galaxies.) The 
FIR-HCN relation is linear all the way up to Lir — 10 12 ' 5 L Q , but the FIR/CO ratio is systematically 
higher for high-luminosity galaxies than for those of normal luminosity. 

One obvious difference between the FIR-CO and FIR-HCN relationships is that the range of CO 
luminosity is only 2 orders of magnitude, whereas HCN luminosity covers 3 orders of magnitude, almost 
the same as the IR luminosity. For LIGs/ULIGs the CO is systematically weak compared with the IR. 
Measurements of CO luminosity of a large sample of ULIGs (Solomon et al. 1997; most of which are not 
in Fig. lb since there are no corresponding HCN observations available) show that they all lie above the 
fitted line in Figure lb. 

Another way to demonstrate these points is to compare the luminosity ratios of Lir/Lrcn and 
Lm/Lco with Lir (Fig. 2) . Lir/Lrcn appears to be nearly independent of Lir, whereas the Lir/Lco 
ratio increases substantially with increasing Lir as established by previous work (e.g., Young et al. 1989; 
Sanders, Scoville, & Soifer 1991). These results are also revealed in Table 2, which summarizes the average 
IR, CO, and HCN luminosities among ULIGs, LIGs, normal spiral galaxies, and all galaxies of the entire 
HCN sample. 

In summary, the IR-HCN correlation is linear and is extremely tight over 3 orders of magnitude 
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in luminosity, when compared with the nonlinear IR-CO correlation. While the high luminosity of 
LIGs/ULIGs requires an elevated SFE of the total molecular gas indicated by Lir/Lco, the SFR per unit 
of dense molecular gas — the SFE of the dense molecular gas indicated by (£ir/£hcn) — is almost constant 
and independent of the IR luminosity or total SFR. The tight relationship between L m and Lrcn will be 
further illustrated by comparing it with other correlations in the following two sections (§3.2 & §3.3). 

3.2. Dense Molecular Gas Fraction (Lhcn / Leo) an d Star Formation Efficiency (£ir/£co) 

Figure 3 shows a significant correlation between HCN and CO luminosities with a correlation coefficient 
R=0.92 (R 2 = 0.85). In all figures, except where Lir is itself explicitly plotted, we distinguish LIGs/ULIGs 
(iiR > 10 11 L Q , filled circles) from the less luminous normal spiral galaxies (open circles) so that any 
systematic difference between these two subsamples can easily be seen. This suggests that the more gas-rich 
galaxies tend to have more dense molecular gas (and vice versa), and thus are more luminous (Fig. 1). The 
correlation is much better for normal galaxies (open circles) than for LIGs/ULIGs (filled circles). For normal 
spiral galaxies there is a tight relationship between the HCN and CO luminosities with a slope of 1.0, shown 
as the fit line in Figure 3. There is also a very small dispersion about the fit with <r(logLncN) = 0.14. This 
good linear HCN-CO correlation is the underlying reason that the FIR-CO correlation is linear for normal 
galaxies and CO luminosity is reasonably good at predicting their SFR (FIR). 

All but three of the 31 LIGs have Lrcn above the linear fit (Fig. 3) determined from normal spirals, 
clearly showing that excess HCN emission compared with CO emission is a characteristic of LIGs/ULIGs. 
An orthogonal fit to all galaxies gives a nonlinear relation: 

logL W K = l.mogLco - 4.79, or ^ = 0.1( - f CO ^ 2 ) 38 . (2) 

Leo 10 iu Kkms x pc 

It is the high ratio of HCN/CO for LIGs/ULIGs that causes the breakdown in the IR-CO relation. 

A high ratio of Lhcn/£co is a distinguishing feature of LIGs/ULIGs. As can be seen from Table 2, the 
average ratio systematically increases as IR luminosity increases. Figure 4 shows the IR luminosity of all 65 
galaxies in the sample as a function of the ratio Lhcn/^co- It is immediately apparent that LIGs/ULIGs 
(Lir > 1O 11 L0) have systematically higher HCN/CO ratios than normal galaxies. However, for normal 
spirals only, although there is some scatter with £hcn/^co = 0.02-0.06, there is no correlation between IR 
and Lhcn/^co All of the most luminous galaxies, i.e., ULIGs, have high £ H cn/£co > 0.09. The most 
striking result is that all ( 27/27) galaxies with Lhcn / -^co > 0.06 are luminous (or ultraluminous). The 
2 most luminous galaxies in our sample Mrk 231 and IRAS 17208-0014 have HCN/CO = 0.23 and 0.25, 
respectively. The molecular ratio Lhcn/^co appears to be 100% successful as an indicator of galaxies with 
infrared starbursts. In § 4 we further discuss the contributions of the disk and nuclear regions to these 
global values of Lhcn and Leo as well as the Lhcn/^co ratio. 

While all galaxies with a high HCN/CO ratio are IR-luminous, it is not true that all LIGs have a high 
dense gas fraction, £hcn/£co- There are a group of seven LIGs in our sample that are IR luminous and 
have high HCN luminosities but normal ratios of HCN/CO. These are gas-rich (CO-luminous) galaxies 
with normal SFE. They are the filled circles on the normal galaxy fit line in Figure 3 and in the upper left 
of Figure 4. Their high Lir is simply the result of a tremendous amount of available molecular gas and a 
normal fraction of dense molecular gas. In this sense they are not starbursts globally. They are globally 
using their molecular gas at a normal SFR. Some of the notable examples in our sample are Mrk 1027, 
NGC 1144 (Gao et al. 1997), NGC 6701, and even Arp 55 (Table 1). Luminous prestarburst galaxies like 
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Arp 302 (Gao 1996; Lo, Gao, & Gruendl 1997), and even early stage galaxy mergers in the initial starburst 
phase with extended CO emission (e.g., NGC 6670, Wang et al. 2001), belong to this category as well. 

Figure 5a shows a correlation between Lir/Lco and Lhcn/Lco with a correlation coefficient R = 0.74 
(R 2 = 0.55, there is little effect on the fit whether those HCN limits are used or not). This suggests that 
the star formation efficiency, SFE (Lir/Lco), depends on the fraction of molecular gas in a dense phase 
indicated by (Lhcn/Lco)- This correlation demonstrates the direct connection between the IR and HCN 
luminosities (as shown in Fig. la). Both Lin and Lhcn have been normalized by Lqo to show the 
physical relationship between Lir and Lhcn after removing the dependence upon distance, galaxy size, 
and other possible selection effects. Of course, the range of the ratio (Lhcn/Lco) is less than that for 
£hcn and the correlation is not as strong as in Figure la, although the dispersion about the fit is almost 
the same. This normalization is dust extinction-free, unlike the use of, for example, the blue luminosity in 
the normalization. 

Similarly, we can show the correlation between Lir and Lco divided by Lhcn for normalization. 
Surprisingly, the strong correlation observed in Figure lb has completely disappeared, and there is no 
correlation left at all (R 2 = 0.01; Fig. 5b). This certainly reflects the tight correlation between IR and 
HCN and the fact that Lir/Lrcn is almost independent of galaxy luminosity (Fig. 2a). The absence of a 
correlation in Figure 5b also implies that the apparent well-known correlation between Lir and Leo may 
lack a true physical basis, at least as compared with that between Lir and Lhcn- The IR-CO correlation 
may be due to a combination of the tighter IR-HCN and HCN-CO correlations. Nevertheless, Lir/Lhcn 
is not completely independent of the CO luminosity (see §3.3 and the appendix) as its inverse Lhcn/Lir 
is weakly correlated with Lco /Lir with a correlation coefficient R = 0.56 (R 2 — 0.31, still a meaningful 
correlation as a consequence of the tight HCN-CO correlation normalized by Lir). 

3.3. The Model Parameter Fits 

Table 3 in Appendix B summarizes almost all the parameter fits in the IR, CO, and HCN data 
sets including the warm dust temperature deduced from the 60 and 100 fim fluxes. Here we discuss the 
three-parameter fits involving the IR, HCN, and CO luminosities and the warm dust temperature Td us t 
dependence of the luminosity ratios. Further details of various other correlation fits are also discussed in 
the Appendix. These fits to the data demonstrate that the two important independent parameters are the 
IR and HCN luminosities. 



3.3.1. Multi-parameter Fits 
The IR luminosity from a model fit to the HCN and CO [the IR(Lhcn, Leo) model] yields 

Zo 3 Lir(Lhcn, L co ) = 2.28 + 0.88?o 3 i H CN + 0.16logL co . (3) 

The weak dependence on Lco shows that the IR luminosity is determined principally from the HCN 
luminosity. The dispersion about the fit is almost the same as the two parameter fit of Lir-Lrcn- The t 
ratio for Lco (0.98) is very small as compared with that for Lhcn (6.95), demonstrating that Lco is 
almost random after fitting Lhcn for Lir. 
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An IR and CO luminosity model fit for Lrcn [the HCN(Lir, Lqo) model] gives 

logL HC N{L m , Leo) = -2.62 + 0.56logL m + 0A9logL GO - (4) 

Here, Lir and Leo seem to be almost equally important. This relation produces a tighter fit than either 
of the two-parameter fits HCN-CO (Fig. 3) or HCN-IR (Fig. la). 

Therefore, although the HCN luminosity depends on both the CO and IR luminosity, the IR is predicted 
almost entirely by the HCN with almost no effect from the CO luminosity. The reason for this surprising 
result is probably that, of the three separate two-parameter (IR-HCN, HCN-CO, and IR-CO) fits, the 
weakest and the one that has the largest dispersion of the data from the fit is the IR-CO correlation. 

The corresponding fit for the CO luminosity is 

logL C o(Lm, £hcn) = 3.60 + 0.55fogL H CN + 0.12logL m . (5) 

It is clear that the CO is more closely related to the HCN than to the IR luminosity. In other words, the 
dependence of the CO luminosity on the IR is much weaker than on the HCN luminosity. This result is 
probably due to the very tight correlation between CO and HCN for normal galaxies (Fig. 3). 

3.3.2. Warm Dust Temperature Dependence 

Here we examine the influence of the warm dust temperature on the luminosity ratios. It is well known 
that the warm dust temperature is an important parameter in the IR-CO correlation (e.g., Young et al. 
1986; Young & scovillc 1991). Calculating the warm dust temperature from the 60/100 ^m flux ratio, we 
find that the Lir/Lco luminosity ratio depends strongly on Tdust (see figures in the appendix) 

Lm/Lco = 61(T dust /35K) 5 - 7 , (6) 

for optically thin dust with emissivity ~ v@ (j3 — 1.5). 

As we can see from Figure 2a, the ratio Lir/Lhcn is almost independent of IR luminosity. If we fit 
this ratio to the dust temperature (poor correlation), we find a weak correlation with 

W^hcn ~ dSOiT^t/SSK) 1 - 8 . (7) 

The IR-HCN relation changes only by a factor of 2 on average across the entire temperature range of 
the sample galaxies (T = 29 to 46K), while the IR-CO ratio changes by a factor of 13. The HCN/CO ratio 
changes with the dust temperature between these two extremes since (using eqs. [6] and [7]) 

Lhcn/Lco = (Lm/ L co ) / (L m / L H cn) ~ T dust . (8) 
Indeed, a direct fit to the data (still a reasonable correlation with R 2 <~ 0.5) gives Lhcn/Lco ~ T^ st . 

3.4. Summary of the Results 

Our results presented in Figures 1-5 show that dense molecular gas mass indicated by HCN luminosity 
is a much better predictor of infrared luminosity and star formation than total H2 content indicated by CO 



-10- 



luminosity. Wc identify two types of LIGs according to their total molecular gas content. Both types have 
high total Lhcn in addition to high Lir. Most LIGs have a high dense molecular gas fraction, but a second 
group has a normal ratio of HCN/CO yet a very high CO luminosity. Are both types of LIGs starbursts? 
It may well be that LIGs with a normal Lrcn / Leo are not necessarily genuine starbursts since the infrared 
luminosity can be produced by a large amount of molecular gas forming stars at a normal rate (efficiency). 
Star formation efficiency can remain almost unchanged, whether galaxies are molecular gas-rich (more 
luminous) or relatively gas-poor (less luminous, cf. Young 1999). However, SFE increases dramatically 
from dense molecular gas-poor galaxies to dense molecular gas-rich galaxies. All galaxies with a dense gas 
mass fraction Lhcn/Lco > 0.06 are LIGs or ULIGs in our sample. In essence, HCN traces the molecular 
gas at high density and at high warm dust temperature that is tightly linked to the active star formation. 

4. DISCUSSION 

4.1. HCN — Tracer of Active Star- forming Dense Molecular Gas 

All the parameter fits essentially show the same thing. CO luminosity by itself leads to a rough 
prediction for IR luminosity that breaks down for luminous infrared galaxies (LIGs), especially for 
ultraluminous infrared galaxies (ULIGs), whereas HCN luminosity is much better at predicting the IR 
luminosity for all galaxies including ULIGs. Therefore, the star formation rate (SFR) indicated by Lir in 
galaxies depends on the amount of dense molecular gas traced by HCN, not the total molecular gas content 
measured by CO. In particular, the IR-CO correlation may not have a solid physical basis as it can be 
readily related to the stronger and perhaps more physical IR-HCN and HCN-CO correlations, which may 
be the origin of the correlation between IR and CO. This is reminiscent of the poor IR-HI correlation as 
compared with the better IR-CO correlation that became apparent two decades ago, when systematic CO 
observations of significant numbers of galaxies became available. 

The HCN radiation is associated with the warm dust traced by the FIR radiation, whereas the 
total molecular gas traced by CO originates in diverse dust components at different temperatures. The 
temperature dependence of Lir (or L e0flm or Lioo^m) is straight-forward to understand since it is just the 
Planck law plus the dust emissivity. In this perspective, it is easy to see why the Lir-Lco relation has a 
strong dependence upon the dust temperature term since Leo at most is proportional to the first power of 
the temperature, while Lir (or Zaoo^m); depending on the emissivity law, is proportional to at least the fifth 
power of the dust temperature. Therefore, many correlations that involve Td us t and Leo [see Appendix, e.g., 
in one case as represented in the CO(L WOlim , Tdust) model, logLco = 5.07 + OSSlogLioo^m ~ SMogT^ust] 
can be easily explained. The more complex question is why this Td us t dependence almost entirely goes 
away for the L H cn, e.g., the HCN(Li 00M m, T dust ) model, ZogL H CN = -1.12 + 0.95logLi 00 ^ m - 0.2logT dusU 
and similarly others involving Lhcn (Table 3). The simple answer to this is that the higher molecular 
gas density produces more active star formation, which raises Td us t owing to heating of the newborn stars. 
High molecular gas density, strong HCN emission, and warm dust temperature go together. HCN traces 
the active star-forming molecular gas where both the molecular gas density and dust temperature arc high. 

We have already shown in Paper I that HCN emission in galaxies is primarily due to the collisional 
excitation by high-density molecular hydrogen, not radiative excitation through the mid-IR pumping (see 
also Stutzki et al. 1988; Paglione et al. 1997). Even though the mid-IR pumping is not a significant source 
to excite the rotational transition of HCN emission, there are still some other possibilities that may help 
excite HCN, e.g., collisions with electrons (Aalto et al. 1995), the possibly enhanced HCN abundance and 
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shock excitation, owing to excess supernovae occurred in starburst galaxies and LIGs/ULIGs. However, 
all these, even collectively, are only a secondary effect on a global scale in galaxies, although significant 
contribution in a particularly favorable environment in some small localized regions cannot be excluded. 

In any case, the physical explanation for the tight correlation between the HCN and IR is star 
formation in dense molecular gas. The active high-mass star-forming sites are the cores of GMCs, where the 
molecular gas is warmer and denser than in the GMC envelopes, where most of the CO emission originates. 
Currently, detailed statistical study examining the relationships among FIR, HCN, and CO on the scale of 
GMCs cores is not yet available. There are extensive observations (Plume et al. 1997) of another dense gas 
tracer, CS emission, of high-mass star formation cores in the Milky Way. All of these regions have an H 2 
density more than sufficient to produce strong HCN emission. 



4.2. Dense Molecular Gas and Star Formation Rates 

For an initial mass function (IMF), typically taken to be the Salpeter IMF, LiR,/M(gas) can be 
interpreted as a measure of star formation efficiency (SFE), i.e., SFR per unit gas mass. This is because the 
SFR is related to Lir by 

M SFR « 2 x l(T lo (Li R /L ) M Q yr-\ (9) 

assuming that the observed FIR emission is produced primarily from dust heating by O, B, and A stars 
(e.g., Scoville & Young 1983; cf. Gallagher & Hunter 1987; Kennicutt 1998b). Although Lir correlates 
with £co, the correlation is nonlinear, with a higher Lir/Lco ratio for higher Lir (Fig. 2b). On the other 
hand, Lir linearly correlates with £hcn, implying an almost constant SFR per unit of dense molecular 
gas mass for all galaxies. 

The HCN luminosity can be related to the mass of dense gas, Mdonsc = o:hcnLhcn, ^ we assume the 
emission originates in the gravitationally bound cloud cores (see Paper I) . For a volume-averaged core density 
n(H 2 ) ~ 3 x 10 4 cuT 3 and brightness temperature T b = 35 K, « H cn = 2.1^n(H 2 )/T fc = 1OM /Kkms _1 pc 2 . 
Substituting in equation (1) gives a luminosity to dense gas mass ratio 

L m /M dense = 9O(a HC N/lO)L Q /M (10) 

for all galaxies, although the mean is actually slightly higher (<~ 120 L Q /M ) for the most luminous galaxies 
(see Table 2). 

Combining equation (1) with equation (9), the SFR in terms of the HCN luminosity is 

M SFR w 1.8 x lO-^LncN/Kkms-ipc 2 ) M^r" 1 . (11) 
In terms of the dense gas mass, the star formation rate becomes 

M SFR «1.8(|^)(^)M j / r- 1 . 

1O 8 M CkHCN 

Since this is a linear relation, the HCN emission is a direct tracer of the SFR in all galaxies 
characteristic depletion time (half-life) is 

Ti/2 = 0.5M dense /MsFR « 2.7a KC NMyrs. (13) 

Although we adopted ckhcn ~ 1OM /Kkms _1 pc 2 for normal spirals (Paper I), this conversion factor might 
be smaller for ULIGs as the T(, can be quite high (Downes & Solomon 1998), leading to shorter dense gas 



(12) 

The dense gas 
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depletion time-scales. HCN observations could potentially become one of the best SFR tracers in galaxies 
in the nearby and distant universe given the high sensitivity and the high spatial resolution achievable at 
millimeter wavelengths with the next generation of the millimeter telescopes. 

There also appears to be some correlation between HCN emission and HCO + , CS, and other tracers of 
star formation, e.g., [C II] line emission and the FIR and 20cm continuum emission (Nguyen-Q-Rieu et al. 
1992), although their sample is very limited. However, it has become clear that [C II] is underluminous in 
ULIGs (e.g., Luhman et al. 2003) and is not a consistent star formation indicator. In addition, in equation 
(9), it was assumed that most of the Lir, originates from star formation with little contribution from active 
galactic nuclei (AGNs) and/or from the general interstellar radiation field. It will be interesting to examine 
correlations between HCN and other indicators of star formation in a large sample of galaxies to assess 
which best indicates the SFR. 

The tight correlation between the IR and HCN emission also implies that the dominant IR emission 
originates from the HCN emission region, especially in LIGs/ULIGs with concentrated molecular gas 
distribution. We know little of the size scales of the FIR emission regions in LIGs/ULIGs. The dominant 
contribution of the radio continuum and mid-IR emission in most LIGs appears to be from the inner regions 
(Condon et al. 1991; Telesco, Dressel, & Wolstcncroft 1993; Hwang et al. 1999; Xu ct al. 2000; Soifcr et al. 
2001). CO emission is usually concentrated in the inner regions, typically within a kpc of the center for 
ULIGs/LIGs (e.g., Scoville, Yun, & Bryant 1997; Downes & Solomon 1998; Sakamoto et al. 1999; Bryant 
& Scoville 1999; cf. Gao et al. 1999). HCN emission originates from the dense cores of the CO emitting 
regions, the sites of star formation, and the source of the FIR emission. Thus, we may predict the size 
scales and location of the FIR emission by determining the HCN source sizes from the HCN mapping. 



4.3. The Global Star Formation Law 

The IR-HCN linear correlation is valid over 3 orders of magnitude from low IR luminosity to the most 
luminous galaxies in the local universe. The direct consequence of the linear IR-HCN correlation is that the 
star formation law in terms of dense molecular gas content has a power- law index of 1.0. The global SFR is 
linearly proportional to the mass of the dense molecular gas (cq. [11] and Fig. 6). Parametrization in terms 
of observable mean surface densities of the dense molecular gas and the SFR will not change the slope of 1 
in the IR-HCN correlation (SFR-Md cn sc correlation), or the linear power index of the star formation law, 
as both quantities are simply normalized by the same galaxy disk area. Our finding of an SFR proportional 
to the first power of the dense gas mass is different from the widely used star formation law with a slope of 
1.4, derived by Kennicutt (1998a) for the disk-averaged SFR as a function of the total (HI and H2) or just 
molecular gas surface density traced by CO emission. As we show below, this law is not valid for normal 
spiral galaxies and results only by combining normal galaxies with starburst galaxies and ULIGs. 



4-3.1. Normal Spiral Galaxies 

The IR-CO correlation (SFR-Mh 2 correlation) is essentially linear up to luminosity L\-r = 10 11 
L Q (ULIGs and most LIGs excluded, see Fig. lb). This seems to also be true in terms of the mean surface 
densities of the SFR and molecular gas mass for the nearest galaxies with spatially resolved observations 
(e.g., Wong & Blitz 2002; Rownd & Young 1999). The linear IR-CO correlation at low to moderate IR 
luminosity is expected since we find that the HCN-CO correlation is extremely tight and linear for normal 
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spiral galaxies (Fig. 3). Thus, the linear form of the global star formation law in terms of total molecular 
gas density as traced by CO for normal galaxies is due to the constant dense gas mass fraction indicated 
by the HCN/CO ratio (discussed in the next section). For normal star-forming spirals, the star formation 
law is linear in terms of both the total molecular gas and the dense molecular gas. Then how did Kennicutt 
(1998a) obtain a slope of 1.4 in the fitting of the star formation law? In our sample, a fit for the normal 
galaxies in the IR-CO correlation gives a slope of 1.0 (Fig. lb), but this is not the case in Kennicutt's 
normal galaxy sample, where there is a poor correlation between the SFR and the gas surface density. 
Thus, no reasonable slope can be derived from his normal galaxy sample alone. 

4-3.2. All Galaxies: Normal Spiral, Starburst, and Ultralurainous Galaxies 

A direct orthogonal regression fit of the IR-CO correlation for all galaxies in our sample (Fig. lb) leads 
to a slope of 1.44. But the best-fitting least-squares slope (errors in Lir, only) is 1.27 (1.25 if galaxies with 
HCN limits are excluded, see §3.1). These fit slopes of the IR-CO correlation are almost identical to the 
star formation power-law index in Kennicutt's (1998a) 36 circumnuclear starbursts and ULIG sample. It is 
obvious from Figure lb that only galaxies with Lir> 10 11 L Q (ULIGs and most LIGs) lie above the fixed 
line of slope 1. This combination of normal and very luminous galaxies leads to a fit with a power-law index 
of 1.4. Therefore, this slope is not a universal slope at all as it changes according to the sample selection. 
The 1.4 slope of the composite Schmidt law (Kennicutt's Fig. 6, 1998a) is determined almost entirely from 
the starburst sample. The circumnuclear starbursts have some of the characteristics of ULIGs/LIGs. In 
particular, they must have a high dense gas fraction indicated by HCN/CO ratios (see section §4.4). 

Indeed, when we add more LIGs/ULIGs into the sample for the IR-CO correlation, the slope becomes 
steeper. In Figure 7, we present an SFR-M(i? 2 ) (IR-CO) correlation diagram with an additional 40 
galaxies, mostly ULIGs, with CO data from the literature (Solomon et al. 1997; Gao & Solomon 1999; 
Mirabel et al. 1990; Sanders et al. 1991), in addition to our HCN sample shown in Figure lb. The 
least-squares fit gives a steep slope of 1.53, and the orthogonal regression fit leads to a much steeper slope 
of 1.73. In this sample luminous galaxies and normal galaxies have equal weight. It is clear that the ULIGs 
steepen the slope of the sample. There also appears to be a trend that some normal spirals with the lowest 
Ssfr and Sh 2 m Kennicutt's sample tend to lie below the 1.4 power fit line. Adding more extreme galaxies, 
both luminous ULIGs and low-luminosity spirals, tends to steepen the slope further toward 2. Therefore, it 
is difficult to derive a unique 1.4 power law based upon the total molecular gas or the total gas content. 

The star formation rate in a galaxy depends linearly on the dense molecular gas content as traced by 
HCN, regardless of the galaxy luminosity or the presence of a "starburst" , and not the total molecular 
gas and/or atomic gas traced by CO and/or HI observations, respectively. Since dense molecular cloud 
cores are the sites of high- mass star formation, it is the physical properties, location, and mass of these 
cores that set the star formation rate. A detailed star formation law can be determined from observations 
directly probing the Milky Way cloud cores, particularly in the Milky Way molecular ring with spatially 
resolved measurements. The molecular tracers that best quantitatively indicate the presence of a starburst 
are primarily abundant molecules with high dipole moments such as HCN and CS requiring high molecular 
hydrogen density for excitation. The molecular property that best characterizes the star formation rate of 
a galaxy is the mass of dense gas. The gas density traced by HCN emission is apparently at or near the 
threshold for rapid star formation. 
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4.4. HCN/CO Ratio A Better Indicator of Starbursts 

Although HCN luminosity is a better indicator of star formation than CO in galaxies, it is very useful to 
compare the HCN with CO to obtain the HCN/CO ratio. The HCN/CO ratio is an indicator of the fraction 
of dense molecular gas available for vigorous star formation and gauges the globally averaged molecular gas 
density. The HCN/CO ratio is also a very successful predictor of starbursts (Fig. 4) and directly correlates 
with the SFE (Lir/Lco) (Fig. 5a). The SFE increases dramatically from dense molecular gas-poor galaxies 
to dense molecular gas-rich galaxies. The global luminosity ratio Lhcn/^co =<SBR> differs among 
galaxies of different luminosities (see Tables 1 & 2) with an average ratio of 1/6 for ULIGs and only 1/25 
for normal spirals. All galaxies in our sample with a global dense gas mass fraction Lhcn/^co > 0.06 are 
LIGs or ULIGs (Fig. 4). We note that Curran et al. (2000) also find an average global ratio of 1/6 in 
luminous IR Seyfert galaxies, which they attribute to star formation. 

The HCN/CO surface brightness ratio (SBR) is potentially a better and more practical indicator than 
the IR/CO ratio (the standard SFE diagnostic) of the starburst strength. Using the ratio IR/CO as a 
diagnostic, the IR emission is presumed to entirely originate from star formation. Other possible sources 
such as AGNs and the general interstellar radiation field are assumed to be negligible. The HCN/CO ratio 
instead is directly related to the molecular gas properties, particularly the local molecular gas density, 
which is tied to star formation. Moreover, the projection and confusion of different velocities along the line 
of sight is inevitably present in IR maps, whereas the different velocity components can be distinguished 
from the kinematics obtained in the detailed CO and HCN maps. In addition, the low spatial resolution 
available in the FIR, even with the Spitzer Space Telescope and the upcoming SOFIA, and other future 
far-IR space missions, is incompatiablc with the high resolution available from millimeter interferometers. 
Therefore, it is important to map the HCN/CO ratio in galaxies in order to fully explore the star formation 
properties and SFE. 

As we show elsewhere from HCN maps of nearby galaxies (Gao 1996; Y. Gao & P.M. Solomon 2004, 
in preparation), the ratio of /hcn/^cOj i-e., the surface brightness ratio (SBR), in most cases, is high in 
the centers of spiral galaxies and typically drops off at large galactic radii. A significant fraction of dense 
molecular gas is (traced by HCN emission) distributed in the inner disks of galaxies outside the nuclear 
or inner ring starburst regions and can be detected to radii as large as a few kpc, perhaps to diameters 
of ~ -D25/4. We find the highest fraction of dense molecular gas, indicated by the SBR (as measured by 
single-dish telescopes with beam sizes <~ 1 kpc), 

SBR (cores) = Ihcn/^co (cores) w 0.1, (14) 

nearly comparable to those observed globally from ULIGs, in the centers of most normal spiral galaxies 
observed (usually the nuclear starburst cores). We attribute this to the presence of a starburst. Heifer & 
Blitz (1993, 1997a) found similarly high ratios in the centers of normal galaxies, which they relate to the 
high ambient pressure, but they do not correlate their data with infrared luminosity or star formation rates. 

We also find that the SBR ratio generally falls off in the disks at larger radii (e.g., ^,3 kpc), to a very 
low SBR ~ 0.015-0.03. This low SBR (disks) £0.03 is the same as that found in the Milky Way's disk 
on average, and over the full extent of nearby GMCs in the Milky Way (Heifer & Blitz 1997b), as well as 
outside the central regions in several normal spiral galaxies mapped with an interferometer (Heifer & Blitz 
1997a). 

Although the global luminosity ratio Lhcn/^co =<SBR> differs dramatically among galaxies of 
different luminosities (see Tables 1 & 2), the difference between the SBR in central beam measurements 
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of galaxies is related to the different telescope beam diameters (e.g., Heifer & Blitz 1993), different source 
sizes (of CO and HCN), and different size scales of the starburst cores in the centers of individual galaxies. 
The CO/HCN intensity ratio (the inverse of SBR) changes from 20 - 80 for normal spiral galaxies to 4 - 
10 for ULIGs. However, if only central beam measurements are used, the ratio seems quite uniform (Aalto 
ct al. 1995). The central beam measurement is not representative of the global measurment (of the entire 
galaxy) and may not even be an accurate measurement of the central region since it is telescope dependent. 
This explains why Aalto ct al. (1995) did not find the strong Lir-Lhcn correlation. 

The distribution of the HCN emission, particularly HCN/CO ratio maps, can be used directly to 
locate starburst sites in nearby galaxies. High-resolution maps of HCN and HCN/CO have been obtained 
in the central regions of a few nearby normal galaxies, some with clear central starbursts. Paglione et al. 
(1995) mapped the central bar in the starburst nucleus of NGC 253 and found a peak HCN/CO=0.2 
falling off to 0.04 at a galactocentric radius of 200 pc. The peak SBR in this starburst is similar to 
the global value in the most luminous galaxies in our sample (Fig. 4). Both the strongest HCN surface 
brightness and the strongest SBR occurred at the location of strong, extended nonthermal radio continuum 
and thermal free-free emission associated with the starburst. Using our Figure 5a, this high SBR 
(0.15-0.2) indicates that this small central starburst has an Lir/£co ratio comparable to that of ULIGs, 
Lyr/ Lco~ 2OOL0/Kkms~ 1 pc 2 . Unfortunately, there are no high-resolution IR observations available to test 
this prediction. Paglione et al. (1995) estimated the Lir from the radio continuum to be £ir^ 1 x 10 9 Lq. 
The Lqo estimated from the HCN and HCN/CO ratios is approximately 4 x 10 7 Kkms -1 pc 2 yielding a 
rough Lir/Lco ~ 25OL /Kkms~ 1 pc 2 , in agreement with the prediction of Figure 5a. 

High HCN/CO ratios (0.1-0.2) are also found in the central £l kpc regions of several other galaxies 
(Downes et al. 1992; Heifer & Blitz 1997a; Reynaud & Downes 1997; Kohno et al. 1996; Kohno, Kawabe, 
& Vila-Vilaro 1999). Kohno et al. (1999) found that most of the HCN emission, as well as an enhanced 
HCN/CO ratio, is associated with the circumnuclear star-forming ring in NGC 6951. No significant 
enhancement of the HCN/CO ratio is observed at the CO peaks. In IC 342 with a resolution of 60 pc 4 
(Downes et al. 1992), however, only three of five HCN clouds seem to be actively forming into stars based 
on the presence of free-free emission. However, the observations of the free-free emission are still of limited 
sensitivities, either at 2cm and 6cm (Turner & Ho 1983; Turner et al. 1993) or at 3 mm (Meier & Turner 
2001). And there are no high-resolution FIR observations available to clearly indicate the star formation 
activities and star formation rates at this small scale. Judging from all these observations and the Ha 
emission (Turner & Hurt 1992), there is probably star formation in all the HCN clouds (D. Downes 2003, 
private communication). Are some of these HCN clouds precursors to a starburst? Although the HCN/CO 
ratios are high (~ 0.14) for these HCN clouds, the average of a several hundred parsec central region has 
only an HCN/CO ratio of ^,0.05 in IC 342, clearly not a strong starburst. 

A recent summary of the HCN/CO ratios obtained in six central nuclear starbursts (Shibatsuka et al. 
2003) lists a typical ratio HCN/CO«0.1-0.25. Sorai et al. (2002) also find the highest HCN/CO ratio of 
^0.1 in central regions of a few nearby galaxies. This is the same as the global ratio we find in many LIGs 
and all ULIGs that have HCN luminosities several hundred times greater than the central starbursts of 
normal galaxies. 

The only LIGs that have been well imaged in HCN are the AGN/starburst hybrid galaxy NGC 1068 
(Tacconi et al. 1997, 1994; Heifer & Blitz 1995; Jackson et al. 1993), which has very high HCN/CO ratio 



4 Distancc to IC 342 is adopted as 3.7 Mpc, see Table 1. 
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(^0.3) within ~ 100 pc nuclear region, the merging pair Arp 299 (Aalto et al. 1997; Casoli et al. 1999), and 
the archetypal ULIG Arp 220 (Radford et al. 1991; N. Z. Scoville 2001, private communication). Arp 299 
may harbor an AGN in the eastern nucleus (e.g., Gehrz et al. 1983), similar to NGC 1068. It appears that 
some weak Seyferts also have rather high HCN/CO ratio (^0.2) in the innermost ,$,100 pc nuclear region 
around the AGN (e.g., M51, Kohno et al. 1996). But this contributes little to the average HCN/CO ratios 
in the circumnuclcar (^0.5-1 kpc) starburst rings where most of the molecular gas is located. 

Ultimately, it is the dense molecular gas, rather than the total molecular gas content, that is the raw 
material for active star formation in galaxies. Our global measurements of 65 galaxies show that the dense 
molecular gas fraction indicated by Lhcn/^co is an important measure of the star formation efficiency; the 
mass of dense molecular gas indicated by Lhcn is a very good measure of the star formation rate deduced 
from Lir. The total molecular content indicated by Leo is an unreliable indicator of star formation rate 
particularly in starburst galaxies. The location of starburst regions and better characterization of the 
starbursts in individual galaxies could be better indicated by the local SBR measurements, obtained by the 
high-resolution and high-sensitivity HCN and CO observations. 



4.5. Starburst Origin of the Far-IR Emission 

Our results show that high-mass star formation in dense molecular gas is responsible for the 
infrared luminosity from a wide range of galaxies including normal spirals of moderate IR luminosity 
(5 x 1O 9 Lq^,Lir^,1O 11 L0), luminous infrared galaxies (LIGs, 10 11 L©^,Lir^,10 12 L q ), and ultraluminous 
galaxies (ULIGs, Lir^1O 12 L ). The star formation rate or IR luminosity is proportional to the mass of 
dense molecular gas in all galaxies. 

4- 5.1. Luminous and Ultraluminous Infrared Galaxies 

Evidence is mounting that the dominant energy source in most ULIGs is from the extreme starbursts 
rather than the dust-enshrouded AGNs (e.g., Solomon et al. 1997; Downes & Solomon 1998; Genzel et al. 
1998; Scoville et al. 2000; Soifer et al. 2001). The debate about the energy source in ULIGs, particularly in 
"warm" ULIGs, has been going on for over a decade (e.g., Sanders et al. 1988; Veilleux et al. 1999; Sanders 
1999; Joseph 1999). Although AGNs are present in many ULIGs (e.g., Nagar et al. 2003; Franceschini 
et al. 2003), there is little evidence indicating that the AGN contribution is the dominant source of the FIR 
emission. 

Our results summarized in § 3 provide compelling evidence in favor of a star formation origin for the 
huge infrared luminosity from ultraluminous galaxies. We have shown that the infrared luminosity of all 
molecular gas-rich spiral galaxies including ULIGs is proportional to the dense star-forming molecular gas 
mass traced by LnCN(see Fig. la and Fig. 6). 

One of the main arguments in favor of an AGN as the power source in ULIGs is the anomalously high 
ratio Lib,/ Leo or Lir/M(H2). The IR luminosity and thus the required star formation rate per solar mass 
of molecular gas (as traced by CO emission) is an order of magnitude higher in ultraluminous and most 
luminous galaxies than in normal spiral galaxies, suggesting that an AGN rather than star formation is 
required (Sanders et al. 1988, 1991) in order to produce the high infrared emission. Viewed in terms of the 
dense gas mass the situation is completely different. The ratio Lir/Lhcn or Lin/Mdense is essentially the 
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same in all galaxies including ULIGs. Figure 8 shows that -km/Mdcnsc is virtually independent of galaxy 
luminosity and on average £iR/Md C nsc — 9OL /Mq, about the same as in molecular cloud cores but much 
higher than in GMCs as a whole (Mooney & Solomon 1988). Ultraluminous galaxies simply have a large 
quantity of dense molecular gas and thus produce a prodigious star burst that heats the dust. It is not 
surprising that starbursts of this magnitude are observed in the infrared and never seen in the optical-UV 
part of the spectrum. The young OB stars are imbedded in very massive and dense regions dwarfing 
anything found in Milky Way GMCs and all of their optical-UV radiation is absorbed by dust. Although 
ULIGs are not simply scaled-up versions of normal spirals in terms of their total molecular mass, they are 
scaled-up versions in terms of their dense molecular gas mass, which is exactly what is expected if star 
formation is the power source. 

Ultimately, even for warm ULIGs (e.g., Surace et al. 1998), which might have, to some degree, evolved 
into the phase of the dust-enshrouded QSOs/AGNs (e.g., Sanders et al. 1989; Surace & Sanders 1999; 
Vcilleux et al. 1999; Genzel et al. 1998; Sanders 1999), we may still be able to tell whether starbursts still 
dominate most of the high Lir or not by examining their total dense molecular gas content and the fraction 
of dense molecular gas. For instance, IRAS 05189-2524 and Mrk 231 are warm ULIGs, but they have 
similar LiR/Md en sc and HCN/CO ratio as the other seven out of nine ULIGs in our sample. Clearly, more 
HCN observations are required to judge if warm ULIGs distinguish themselves from other ULIGs. 

4-5.2. Normal Spiral Galaxies 

Although some fraction of Lir originates from outer disks in nearby large spiral galaxies, the dominant 
contribution is still from the centers (e.g., Rice et al. 1988). Better resolution IRAS maps, obtained by an 
improved imaging deconvolution algorithm, tend to show much more centrally concentrated FIR emission 
in the inner disks (Rice 1993). Devereux & Young (1990) have shown that the global FIR luminosity in 
spiral galaxies is consistent with that contributed by warm dust heating from the high-mass OB stars. 
Higher resolution (as compared with that of IRAS ) measurements of the 160 /xm and Ha emission in 
nearby galaxies NGC 6946 and M51 suggest that the FIR luminosity is in quantitative agreement with that 
expected from OB stars throughout the star- forming disks (Devereux & Young 1992, 1993). 

The existence of cold dust components of <J 30 and ^ 20 K in nearby galaxies has been revealed from 
both recent ISOPHOT FIR observations (e.g., Popcscu et al. 2004; Haas et al. 1998; Alton et al. 1998a; 
Davies et al. 1999) and SCUBA submillimeter observations (e.g., Alton et al. 1998b; 2000), respectively. 
The cold dust distribution usually has a larger radial extent (e.g., Alton et al. 1998a) and the dominant 
cold dust component may not be closely associated with the active star-forming inner disks. HCN maps of 
nearby galaxies (e.g., Nguyen-Q-Rieu et al. 1992; Heifer & Blitz 1997a; Y. Gao & P.M. Solomon 2004, in 
preparation) show that the HCN emission region is much more compact than that of CO emission region in 
normal spiral galaxies and/or starburst galaxies and is thus closely related to the warm dust in the inner 
disks. 

The contribution of the general interstellar radiation field to the total FIR emission in galaxies might 
be significant at the Iow-Lir end, where the general infrared interstellar radiation field is comparable to 
the IR radiation from the active star formation. This might be testable on the IR-HCN correlation plot 
(Fig. la) for galaxies of the lowest Lir (^5x 10 9 L q ) when more observations of the lowest Lir sources 
are available to show a statistically significant trend. Given that the tight correlation between these two 
quantities might be used to predict one another, one can check whether Iow-Lir galaxies have a higher 
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Lir/Lhcn than expected from the IR-HCN correlation. 

4.6. High Redshift Galaxies and AGN 

AGNs may become more important than star formation at the very high IR luminosity end, especially 
in warm ULIGs (Vcilluex et al. 2003), in contributing the bulk of the energy output (e.g., Sanders et al. 
1988; Veilluex et al. 1999). For extremely luminous galaxies with Lir~ 10 13 L q , the so-called hyperluminous 
infrared galaxies (HLIGs, Sanders & Mirabel 1996; Rowan-Robinson 2000), the implied SFR (Equation 9) 
would be ~ 2000 M Q yr _1 . Although such super-starburst galaxies likely exist at high redshifts and HLIGs 
could indeed be such super-starburst systems (if the magnification by a gravitational lens and the AGN 
contribution to the IR emission are not important), there are no similar systems in the local universe. 

If the tight correlations obtained from our local HCN sample are applicable to high- 
z galaxies, we can roughly estimate their expected molecular gas properties. For the whole 
sample LiR,/LHCN=900L Q /Kkms _1 pc 2 (Equation 1). But ULIGs have a slightly higher ratio of 
£iR/-kHCN=1200L Q /Kkms~ 1 pc 2 (Table 2). For high-z galaxies with Lir^ 10 13 L q , we thus expect to 
have Lhcn ~ 0.8 x 10 10 Kkms _1 pc 2 if they are the analogs of local ULIGs. This is only a factor of 2-3 
higher than the highest Lhcn observed in local ULIGs. The first high-z HCN detection is from the 
Cloverleaf quasar at z=2.567, helped by the magnification of the gravitational lens, with the intrinsic, 
magnification-corrected L H cn ~ 0.3 x 10 10 Kkms _1 pc 2 (Solomon ct al. 2003). It appears that, if the hot 
dust AGN component can be subtracted from the total IR emission, then even the Cloverleaf quasar seems 
to fit the IR-HCN correlation determined here from the local universe. 

At present, more than 20 CO sources at high z have been detected, all with high CO luminosity 
(e.g., Cox et al. 2002; Carilli et al. 2002; Papadopoulos et al. 2001; Guilloteau et al. 1999; Downes 
et al. 1999; Frayer et al. 1998). We can estimate the expected Lhcn for these extraordinarily large 
Leo sources detected at high z using our HCN-CO correlation (Fig. 3). Cox et al. (2002) report the 
strongest CO emitter detected to date, with Leo ~ 10 11 Kkms _1 pc 2 even after the magnification by the 
gravitational lensing is corrected. Using Equation 2, or Lhcn/Lco = 0.1 x (Lco/10 10 Kkms _1 pc 2 ) - 38 , 
we obtain Lhcn ~ 2.4 x 10 10 Kkms~ 1 pc 2 and Lhcn/Lco=0.24. This Lhcn/Lco ratio is the same as 
the very highest found for local ULIGs. For this largest Leo source, L IR ~ 2.0 x 10 13 L Q is estimated by 
Cox et al. (2002) based on the 250 and 350 GHz measurements (Omont et al. 2001; Isaak et al. 2002). 
And using LiR/LHCN=1200L Q /Kkms~ 1 pc 2 , we expect to have Lhcn = 1-7 x 10 10 Kkms _1 pc 2 . Perhaps 
a better estimate of the expected HCN luminosity can be constrained from both IR and CO luminosities 
by using the multiparameter fit from equation (4). Indeed, we also obtain Lhcn = 1-7 x 10 10 Kkms^" 1 pc 2 . 
If such high Lhcn is eventually detected in this source, then star formation rather than AGNs must be 
responsible for most of the high infrared luminosity. 

We can also roughly estimate the expected upper limit of Lir/Lco f° r HLIGs and/or high-z galaxies, 
if they are powered by star formation and the correlation between Lir/Lco and Lhcn/Lco found in 
Figure 5a is applicable. Although there are very large scatters (~ 0.5 dex, 2a) in the fit, this should be 
good to within a factor of <~ 3. Here we can take Lhcn /Leo ~1 as upper limits. Using the orthogonal 
fit logL m /L C o = L24;o#Lhcn/Lco + 3.24, we obtain Lni/Lco^^OOLQ/Kkms^pc 2 . Therefore, the 
expected maximum Lir from star formation is always less than £1.7 x 10 14 (Lco/10 11 Kkms -1 pc 2 )L Q . 

Although many submillimeter galaxies are HLIGs, few have Lir ~ 10 14 L Q (Chapman ct al. 2003). The 
highest Leo detected among submillimeter galaxies is ~ 0.7 x 10 11 Kkms _1 pc 2 (Ncri ct al. 2003; Greve, 
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Ivison, & Papadopoulos 2003). If any sources indeed have high Lhcn^ 0.7 x 10 n Kknis _1 pc 2 , it appears 
that extreme starbursts from abundant active star-forming dense molecular gas are still possible to power 
such L IR ~ 10 14 L Q sources. This corresponds roughly to the maximum possible SFR of ~ 2 x 10 4 M Q /yr 
(using eq. [9]), as argued by Heckman (2000) based on physical causality for self-gravitating, extremely 
massive galaxy systems of ~ 10 x L*. Observing HCN in addition to CO will help distinguish whether star 
formation rather than AGNs contributes mostly to the highest Lir, observed. 

5. SUMMARY 

HCN luminosity is a tracer of dense molecular gas, n(H 2 )^3 x 10 4 cm~ 3 , associated with star-forming 
molecular cloud cores. Here we briefly summarize the principal results found from an analysis of our HCN 
survey of galaxies: 

1. A tight linear correlation between the IR and HCN luminosities Lir and £hcn in 65 galaxies is 
established with a correlation coefficient R=0.94 (Fig. la). There is also a significant correlation between 
the normalized luminosities Lir/Lco and Lhcn/Lco which confirms the true physical relationship between 
Lir and Lhcn • The IR-HCN linear correlation is valid over 3 orders of magnitude including ultraluminous 
infrared galaxies, the most luminous galaxies in the local universe. The direct consequence of the linear 
IR-HCN correlation is that the star formation law in terms of dense molecular gas content has a power 
law index of 1.0. The global star formation rate is linearly proportional to the mass of dense molecular gas 
in normal spiral galaxies, luminous infrared galaxies and ultraluminous infrared galaxies. This is strong 
evidence in favor of star formation as the power source in ultraluminous infrared galaxies since the star 
formation in these galaxies appears to be normal and expected given their high mass of dense star-forming 
gas. 

2. The star formation rate indicated by Lir depends on the amount of dense molecular gas traced by 
HCN emission, not the total molecular gas traced by CO emission. One of the main arguments in favor of 
an AGN (Sanders et al. 1988) as the power source in ultraluminous infrared galaxies is the anomalously 
high ratio Lir/Lco or £ir/M(H 2 . The IR luminosity and thus the required star formation rate per solar 
mass of molecular gas traced by CO emission, is an order of magnitude higher in ultraluminous infrared 
galaxies and most luminous infrared galaxies than in normal spiral galaxies. This has been interpreted 
as indicating that a dust enshrouded AGN rather than star formation is required to produce the very 
high luminosity (Sanders et al. 1991). Viewed in terms of the dense gas mass the situation is completely 
different. The ratio Lir/£hcn or -^IR/Mdcnsc is the same in all galaxies including ultraluminous infrared 
galaxies. Figure 8 shows that LiR/M^ense is virtually independent of galaxy luminosity and on average 
iiR/Mdcnso ~ 90Lq/M q , about the same as in molecular cloud cores but much higher than in GMCs as a 
whole (Mooney & Solomon 1988). Ultraluminous infrared galaxies simply have a large quantity of dense 
molecular gas and thus produce a prodigious starburst which heats the dust, produces the IR, and blocks 
all or most optical radiation. Although ultraluminous infrared galaxies are not simply scaled up versions 
of normal spirals in terms of their total molecular mass they are scaled up version in terms of their dense 
molecular gas mass, which is exactly what is expected if star formation is the power source. We note 
that our sample includes 9 ultraluminous infrared galaxies (£ir> 0.8 x 10 12 L Q ) and 23 luminous infrared 
galaxies (£ir> 10 11 L©) but only 2 "warm" ultraluminous infrared galaxies (Surace et al. 1998), a subclass 
representing about 25% of all ultraluminous infrared galaxies. Although these 2 warm ultraluminous 
infrared galaxies appear normal in terms of the IR/HCN ratio, further investigation of warm ultraluminous 
infrared galaxies is required. 
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3. The HCN/CO ratio is an indicator of the fraction of dense molecular gas available for vigorous star 
formation and gauges the globally averaged molecular gas density. It is a powerful starburst indicator. 

There is a strong correlation between the HCN and CO luminosities in galaxies although the 
correlation is not linear. The ratio £hcn/£co is constant for normal spirals and increases for luminous 
and ultraluminous IR galaxies. The global ratio Lhcn/^co =<SBR> (see Tables 1 & 2) has an average 
ratio of 1/6 for ultraluminous infrared galaxies and only 1/25 for normal spirals. The HCN/CO ratio is 
a very successful predictor of starbursts (Fig. 4) and directly correlates with the star formation efficiency 
indicator (£ir/ Lqo) (Fig. 5a). The SFE increases dramatically from dense molecular gas-poor galaxies to 
dense molecular gas-rich galaxies. All galaxies in our sample with a high dense gas mass fraction indicated 
by Lucn/Lqo > 0.06 are luminous or ultraluminous infrared galaxies (Fig. 4). 

4. The correlation between Lir and Lqo may be a consequence of the stronger correlations between 
-Lir, and Lhcn and between Lhcn and Lqo- A model two parameter fit for -^ir,(£hcn > Lqo) shows almost 
no dependence on Leo- Much of the CO emission originates from moderate density regions with low to 
moderate dust temperature and little or no active high mass star formation. A critical molecular parameter 
that measures star formation rates in galaxies is the amount of dense molecular gas measured by the HCN 
luminosity. High molecular gas density, strong HCN emission and warm dust heated by the newly formed 
OB stars go together. 

5. A quantitative star formation theory must start with the processes that form dense cloud cores, 
particularly massive cores. It appears from our survey of 65 IR/CO bright galaxies including normal spirals 
and very luminous IR starbursts that once the local density is raised from that of a typical CMC (n(H 2 ) ~ 
a few xlO 2 cm -3 ) to 3 x 10 4 cm~ 3 , star formation including high-mass star formation is efficient and 
progresses rapidly. The characteristic time scale for using half of all the dense gas is about 2 x 10 7 years. 

Many thanks are owed to Judy Young for helpful discussions and suggestions. We appreciate the 
generous support and allocation of observing time from the NRAO 12m, IRAM 30m, and FCRAO 14m 
telescopes. We also thank the anonymous referee for a careful and helpful report. This research has made 
use of the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propulsion Laboratory, 
Caltcch, under contract with the National Aeronautics and Space Administration. 
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This table contains all HCN Survey data of Gao & Solomon (2003, Paper I) and includes a dozen galaxies 
in the literature (almost entirely from Solomon et al. 1992, but M51 and NGC 4945 from Nguycn-Q-Ricu 
ct al. 1992 and Henkcl ct al. 1994, respectively). LIGs with Lir> 1O 11 L are in boldface, and ULIGs with 
LirXL1O 119 L0 are further marked with an asterisk (). 

b The 2a upper limit (<) is listed for NGC 1055. The lower limits (>) are for nearby galaxies, where we 
only detected HCN in the galaxy centers or more extensive mapping is still required. 

c The warm dust temperature Td us t is for optical thin dust with emissivity — = 1.5. 
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1 Galaxies with only limits in HCN luminosity are excluded. 

2 The inverse of average surface brightness ratio <SBR> = £HCN/ico indicating the fraction of 
dense molecular gas in galaxies. 
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Fig. 1. — (a) The correlation between HCN and IR luminosities in 65 galaxies. Some limits in HCN 
luminosities are indicated with arrows, (b) The correlation between Lco and Lir for the same HCN sample. 
The sample is divided into luminous infrared galaxies (LIGs) and ultraluminous infrared galaxies (ULIGs) 
with Lir > 10 11 L Q and less luminous "normal" spiral galaxies. The solid lines are the fits to the less 
luminous galaxies with a slope fixed at unity. A single slope of 1.0 fits HCN data for both low and high IR 
luminosities, but not for the CO data. 
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Fig. 2. — (a) The correlation between IR/HCN and IR is almost of nonexistence; (b) whereas the correlation 
between Lir/Lco and £ir is rather prominent. The different trends are particularly obvious at high IR 
luminosity for LIGs and ULIGs with Lir, > 1O 11 L . 
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Fig. 3. — Strong correlation exists between Lhcn and Lqo indicating that the more molecular gas- rich 
galaxies tend to have larger amount of dense molecular gas as well. Similar to Figure lb, the fit line is for 
less luminous normal galaxies (open circles) with a fixed slope at unity. Apparently, almost all LIGs and 
ULIGs (filled circles, Lir > 10 11 L Q ) lie above the line. 
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Fig. 4. -— Lhcn/Lco stays at fairly small values for normal spiral galaxies, but increases dramatically 
for LIGs and ULIGs (Lin > 10 11 L Q ). All galaxies with Lhcn/-^co~0.06 in the sample are luminous and 
ultraluminous galaxies. 
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Fig. 5. — (a) Correlation between Lhcn/Lco and Lir/Lco revealing the true relationship between the HCN 
and IR since both luminosities are normalized by Lqo- This removes all dependence on distance and galaxy 
size and shows that there is a true physical correlation between the IR and HCN luminosities. The best fit 
has a correlation coefficient of 0.74. (b) No correlation between -Ltr/Lhcn an d Lco/Lhcn (a correlation 
coefficient of 0.1) suggests that the correlation between Lir and Leo observed in Figure lb may not be a 
truly physical relation. The sample is divided into LIGs and ULIGs with Lir > 10 11 L Q (filled circles) and 
less luminous normal spiral galaxies (open circles). 
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Fig. 6. — The global star formation law of dense molecular gas. The star formation rate is linearly 
proportional to the dense molecular gas mass. The solid circles are for LIGs and ULIGs with Lir^1O 11 L0, 
whereas the open circles are for the less luminous normal spiral galaxies. 
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Fig. 7. — The star formation rate vs. total molecular gas mass. The solid circles are for LIGs and ULIGs 
with Lir^1O 11 L , whereas the open circles are for the less luminous normal spiral galaxies. More CO data 
of mostly ULIGs available from literature (stars) are added to our HCN sample (see Fig. lb). The orthogonal 
least-squares fit now has a slope of 1.73 with a correlation coefficient R — 0.89. The line of fixed slope of 1, 
a valid fit for normal spirals, is also shown for comparison. 
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Fig. 8. — The star formation rate (far-IR luminosity) per unit of dense molecular gas mass is essentially 
independent of the total IR luminosity. The average for the entire sample is £iR,/Md C nsc = 9OL /M Q . 
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A. Correlations with the Warm Dust Temperature 

The temperature T dust can be estimated as a function of /eo^m/ '/loOjum and where the emissivity 
is the Planck function of a single warm dust temperature times the frequency to the power of /3, i.e., v$ . 
This is for warm dust Tdust between ~ 25 and 60 K and (3 ~ 1 — 2 (e.g., see the appendix of Lonsdale 
ct al. 1985). Although there is increasing evidence for a substantial amount of cold dust in spiral galaxies 
(e.g., Alton et al. 2000), the warm dust is probably still the most important component, especially in 
luminous and ultraluminous infrared galaxies (LIGs and ULIGs), and often a single dust temperature 
fitting of the far-IR spectral energy distribution (SED) including submillimeter measurements is still a good 
approximation (e.g., Lisenfcld, Isaak, & Hills 2000; Dunne et al. 2000). 

The tightness of the IR-HCN correlation can be further examined by comparing the ratio of Lir/L H cn 
with feovun/ fiootim- Figure 9a shows that Lir/L H cn has only a weak dependence upon / 6 o M m/ '/loo^m or 
Tdust, whereas Figure 9b indicates that Lir/Lco correlates strongly with feo^m/ /ioo/xm (R 2 — 0.72) or Td us t- 
The fits in terms of dust temperature have been given in equations (6) and (7). Unlike Figure 5b, where 
Lir/Lhcn appears to be totally independent of Lco/^hcn, Figure 9a indicates that Lir/Lhcn might still be 
weakly dependent upon Tdust- There is also a meaningful correlation between Lhcn/Lco and f&o^m/ '/ioo/xm 
(R 2 = 0.46). An orthogonal fit can be obtained as Lhcn/Lco = 10~ - 71 x (/eo^m/ /loo^m) 2 ' " 1 " ' 8 with a 
fairly large uncertainty in the slope since the correlation is not tight. 

It is not surprising that Lir/Lco is a strong function of Td us t since the dust radiates thermally as 
^dust to ^dusf depending on the dust emissivity to produce IR emission (e.g., Soifer et al. 1989). Thus, 
Figure 9 is not independent of Figure 2 as Td us t is closely related to Lir. And Leo depends only linearly 
on the temperature, ~ Td ust , if the dust and gas are coupled and the gas is thermalized with the intrinsic 
brightness temperature T& ~ T^ust ( e -S-> tne black-body model of Solomon et al. 1997), and more weakly 
in the power of Td us t dependence if not coupled/thermalized. HCN luminosity seems also to depend, at 
a first glance, only on the first power of Td us t if T& ~ T^ ust , but it is actually not simply proportional 
to T^ ust since the HCN emission does not fill the source area probed by the telescope beam, whereas 
CO has probably an area filling factor close to unity, especially for LIGs/ULIGs. Recent interferometric 
HCN imaging in some nearby galaxies indeed shows much smaller HCN source size than that of CO (e.g., 
Downes et al. 1992; Heifer & Blitz 1997a; Kohno et al. 1996, 1999). It is possible that the HCN filling 
factor depends on some power of Td us t- In addition, the weak correlation between Lir/Lhcn and Td us t 
indicates that there might be other important parameters, probably the molecular gas density, responsible 
for the high Lhcn in LIGs/ULIGs. Thus, Lhcn depends on a much higher power of Td us t, and the ratio of 
Lhcn/Lco should depend fairly strongly on some power of Td us t (eq. [8]) rather than being independent. 

This rather strong Td us t dependence of Lhcn can be understood since Lhcn /Leo also correlates with 
Lir (Fig. 4). Because the Lhcn/Lco ratio is such a good indicator of starburst and most LIGs/ULIGs 
have the highest Lhcn/Lco ratio, warmer LIGs/ULIGs of higher Tdust (higher Lir) should have a higher 
fraction of the dense molecular gas (higher Lhcn /Leo), and thus very likely high average gas density. 
Therefore, Lir somehow correlates with the molecular gas density and Td us t is thus also related to the gas 
density, which should be expected as the gas and dust are somehow coupled. 
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B. Other Correlations and Model Fits 

B.l. Simple Two-parameter Fits 

Although Lhcn and Leo are strongly correlated (Fig. 3), there are significant differences, however, 
when they are compared with either -Lioo^m or Lir (Table Bl). The better correlation between Lhcn and 
Lir (vs Leo and Lir) can be recognized in both R 2 (§3.1) and the R.M.S. deviations of the correlation 
fits (Table 3, 0.24 vs 0.34). The luminosity Lioo/im also correlates better with Lhcn than with Leo as the 
difference in R 2 is also significant (0.89 vs 0.80). Nevertheless, Leo appears to be much better correlated 
with Lioo^m than with Lir (there is a significant difference in the logarithmic R.M.S. deviations: 0.22 vs 
0.34), though the difference in R 2 (0.80 vs 0.77) is small. In comparison, there is little difference in the 
correlations between Lioo^m and Lhcn vs between Lir and Lhcn (i? 2 =0.89 vs 0.87 and logarithmic 
R.M.S. 0.22 vs 0.23). But the HCN correlations are much tighter than the CO correlations, even though 
the Lioo M m~Lco correlation has same R.M.S. deviation. 

The luminosity Lioo^m traces relatively the cooler component of the warm dust emission at Tdust ~ 25- 
50 K. In this regime the 100 [im emission is a reasonable tracer of the warm dust mass Md us t- In comparison, 
the total IR luminosity Lir contains the mixture of the various dust components including the hot dust 
radiated at the mid-IR emission. We here estimate the warm dust mass Md us t from the 100 fim flux density 
and the warm dust temperature Td us t and correlate it with Lhcn and Leo in Table 3 as well. It appears 
that Leo is only slightly better correlated with Md us t than Lhcn with Md us t- In short, the differences in 
the above-mentioned various correlations suggest that HCN is a much better tracer of both the total IR 
and 100 ^m emission than CO, whereas CO only traces better the 100 /im than the total IR emission. The 
dense molecular gas rather than the total molecular gas is more intimately related to both the total IR and 
100 /im. emission. 

B.2. Multi-parameter Fits 

The predicted Lir from the CO and Td us t [the IR(Lco> ?dust) model, logLi^Lco, Tdust! model. fit) = 
—6.34 + 0.97logLco + 5.4ZogTd us t] turns out to be the tightest correlation (R 2 = 0.94) among all. The 
prediction of Lir from CO and Td us t is more accurate than from HCN alone and also slightly better than 
from HCN plus Td us t- This surely suggests the importance of Td us t when the total molecular gas content, 
rather than the dense molecular gas, is concerned in predicting Lir. This also implies that Td us t is an 
important parameter in regulating star formation of the total molecular gas, but not necessarily the dense 
(active star-forming) molecular gas. 

We also use Td us t as well as Lioo^m and L &0 ^m, though they are not independent, in the three-parameter 
model fits as listed in Table 3, and we find the following: 

1. The value of Leo decreases with Td us t for given Lioo^m or Lir. The value of Lhcn is, however, 
linearly proportional to Lioo^m or Lir with no correction or only marginal correction of Td us t- Thus, 
Lhcn decreases little with Tdust for a given Lioo^m or Lir. To put it in a different perspective, Lir is 
only slightly better predicted when Td us t is considered together with Lhcn, whereas Lir is much better 
predicted when both Leo and Td us t are in the fit. 

2. The HCN(Lcoj 7d us t) model implies that HCN correlates with Td us t to the 3.3 power as well as 
linearly (the first power) with Leo (cf. Equation 8). Thus, HCN responds sensitively to Td us t with a 



Table 3. Model Fit Parameters 



Para. L 60M m Li o M m Ltr T du3t L H cn L C o Const. RMS a R 2h 



3 parameter correlations 



LlR, 


0.88 0.16 


2.29 


0.24 


0.87 


Lm 


5.4 0.97 


-6.34 


0.17 


0.94 


Lin 


2.9 0.86 


-0.48 


0.20 


0.92 


Lhcn 


3.3 1.00 
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0.21 


0.89 


Lhcn 


0.56 0.49 
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0.96 (-1-7) 
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simple 2 parameter correlations 
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0.72 


0.34 


0.77 


Lm 


1.02 


2.58 
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0.88 


Lhcn 


1.26 
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0.28 
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Lhcn 


0.91 


-1.93 


0.23 


0.87 


Lhcn 0.96 




-1.52 


0.22 


0.89 


Lco 0.72 




-2.12 


0.22 


0.80 
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0.62 


2.49 


0.30 


0.75 


Lco 


0.67 


3.92 


0.27 


0.85 


M dust 


0.84 


0.41 


0.25 


0.84 


M dust 


1.13 


-3.35 


0.21 


0.87 


a rms = Root Mean Square deviation (lo 


garithmic) between the data (observed) and the predicted from the (correlation) 


model fits. 










b The squared correlation coefficient R? . 










Note. — All quantities are logarithmic in 


the correlation fits. The line luminosities are 


in units of Kkm s 


1 pc 2 and the total 


IR, 60, and 100 ^m luminosities are in units of Lq, while Td us t is in K. The warm dust 


mass Mdust was estimated 


using only 


the 100 (im flux and warm dust temperature. The model fit for Lm from Lhcn and Lco (the IR(LhcN; 


Lco) model), for 



example, is logLm{p re dicted) = 0.88Zo</Lhcn + 0.16io<?Lco +2.3 with the corre lation coefficient R = 0.93 (R? = 0.87) and the 
rms error is (E[logLm(p re dicted) — logLm.(obs.)] 2 / N) 3 = 0.24. Numbers with parentheses indicate only marginal significance. 
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power of ^4.3. The IR(Thcn, Tdust) model gives a 2.9 power correlation with Td us t and a 0.86 power with 
LhcNi again roughly consistent with the high (~ 6) power Tdust dependence (cf. eq. [6]) when the Tj U st 
dependence of Lrcn (eq- [8]) has been considered. In short, the difference in the dependences of the HCN 
and IR upon the warm dust temperature is not dramatic. Thus, the ratio of Tir/Thcn has only a weak 
dependence upon Td us t (eq. [7] and Fig. 9a). 

3. The HCN(i60Aim, Tioo^m) model fit and 00(1/60/™, Tioo^m) 3 parameter correlation model reveal 
that HCN correlates primarily with Tioo/im, and essentially not at all with ieo^m, whereas CO correlates 
with both Leo^ m and Tioo/jm- This again implies the much tighter dependence of CO upon Tdust, and is 
basically reiterating the first point described above since there are only two independent parameters among 

Tdust, Tgo^m, LlOOjum an d Tir. 



B.3. Fits Between the Ratios 

Figures. 4 and 5 already showed some correlations between the parameter ratios. Here we just list the 
weak correlation between Tir/Thcn and feo^ni/ fioopim (i? 2 =0.17, Fig. 9a), as directly obtained from the 
orthogonal fit, Tir/Lhcn = 10 3 1 x (feo^m/ fioa^m) ' 63 with a logarithmic R.M.S. of 0.21. The T dus t power 
and the fit have large scatters owing to the poor correlation. The scatter is about the same or only slightly 
better than the R.M.S. of the observed mean luminosity ratio (in logarithm) of log <T ir /L H cn> f° r the 
entire sample, which is ~ 0.23. The correlation fit for L m / Lqo vs f&Oum/ fwop.m or T dust (Fig. 9b, eq. [6]) is 
significantly tighter and has a logarithmic R.M.S. of 0.17. This small R.M.S. is obviously much better than 
the R.M.S. of the observed mean log <Tir/Lco> of the entire sample, which is ~ 0.37. 
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Fig. 9. — (a) Lir/Lhcn as a function of the far-infrared (FIR) color (the 60/um-to-100/zm flux ratio) which 
indicates the warm dust temperature Td us t- (b) Lir/Lco as a function of the FIR color (or Tdust)- There is 
only weak correlation between Tir/Thcn and Tdust, whereas strong correlation exists between Lir/Lco and 
Tdust- For optically thin dust with emissivity oc v 1 - 5 , Td us t is approximately in a range of 25-50 K for galaxies 
in our sample. The sample is divided into luminous and ultraluminous infrared galaxies (LIGs/ULIGs) with 
Tir > 10 11 L Q (filled circles) and less luminous galaxies (open circles). 
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